Astronomy & Astrophysics manuscript no. Franz 


© ESO 2009 


September 25, 2009 





On the velocity field of sunspot penumbrae 



I. A global view 
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ABSTRACT 

Aims. We investigated the vertical penumbral plasma flow on small spatial scales using data recorded by the spectropolarimeter of the 
solar optical telescope onboard Hinode. 

Methods. We computed maps of apparent Doppler velocities by comparing the spectral position of the Fe I 630.15 nm & Fe I 630.25 
nm lines with the averaged line profiles of the quiet Sun. To visualize the flow pattern in the low photosphere, we used a bisector of 
the wing of the absorption lines. Due to the small heliocentric angle (3° < < 9°) of our data sets, the horizontal component of the 
Evershed flow (EF) does not contribute significantly to the line shift. 

Results. We found that in the quiet Sun (QS), the area showing up-flows is always larger than the one exhibiting down-flows. In the 
penumbra, up-flows dominate only at low velocities |Vdop| < 0.4 km s"', while at larger velocities |Vdopl > 0.6 km s"' down-flows 
prevail. Additionally, the maximal up-flow velocity in penumbrae is smaller, while the maximal down-flow velocity is larger with 
respect to the QS velocities. Furthermore, on a spatial average, the penumbra shows a red-shift, corresponding to a down-flow of more 
than 0.1 km s"'. Up-flows are elongated and appear predominately in the inner penumbra. Strong down-flows with velocities of up 
to 9 km s"' are concentrated at the penumbra-QS boundary. They are magnetized and are rather round in shape. The inner penumbra 
shows an average up-flow, which turns into a mean down-flow in the outer penumbra. The up-flow patches in the inner penumbra 
and the down-flow locations in the outer penumbra could be interpreted as the sources and the sinks of the EF. We did not find any 
indication of roll-type convection within penumbral filaments. 
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1. Introduction 



Although penumbrae of sunspots have been subject to scientific 
investigation for a long time, there remain a lot of open ques- 
tions regarding e.g. their filamentary fine structure, the energy 
transport and the plasma dynamics on scales of less than 1". 

The gappy model (iSpruit & Scharme j2006l) assumes the ex- 
istence of field-free, radially aligned gaps below the t = 1 level 
intruding into a potenti al field above. Dark cores of penumbral 
filaments jScharmer et aL 2002) are explained as a variation of 
' the T = 1 level across the tip of the field-free gap, which is 
similar to the explanation of dark lanes in umbral dots given by 
ISchtissler & V ogler (2006). Even though the gappy model pro- 
vides an explanation not only for the dark cores but also for the 
brightness of the penumbra, it poses a pro blem as it does no t 
explain the magnetized Evershed flow (EF) (iRezaei et al ] l2006l) . 

In the uncombed model jSolanki & Montavonll 19931) . nearly 
horizontal flux tubes are assumed to be embedded in a 
less in clined magnetic background field (ISchlichenmaier et al.l 
ll998 a'!F). Strong gradients of atmospheric parameters are en- 
countered for any line of sight (LOS) penetrating both the 
background atmosphere and the flux tube. It has recently been 
shown that dark cores of penumbral filaments can be explained 
within the framework of this model by taking the hot EF, 
which is embedded in a stratified atmosphere, into account 
dRuiz Co bo & Bellot Rubio 2008). Despite the success of the 
uncombed model in explaining e.g. the EF as well as the ex- 
istence of bright filaments and their migration, the permanent 



stability of magnetic flux tubes in the penumbral photosphere is 
questionable. 

Recently, significant progress regarding simulations of entire 
sunspots has been made. Different 3D MHD c odes were u sed 
to model sunspots in slab ge ometry dH eineman ri et al]|2007[) or 
even the entire spot at once (Re mpel et al.ll2009 ), which repro- 
duces a variety of observed features, e.g. dark cored penumbral 
filaments and their migration. In order to either endorse or dis- 
prove these results and to gain a deeper understanding of the 
convective nature of the energy transport in the inclined mag- 
netic field of the penumbra, it is thus crucial to investigate the 
plasma flow in the low photosphere on scales of less then 1"). 

So far velocities are either inferred f r om in v ersion (e.g. SIR, 
Milne-Eddington) (iBellot Rubio et al.l 120071; IScharmer etal] 
2008) or magnetograms that are constructed in the far wings 
of magnetic absorption lines, thereb y identifying str ongly blue- 
or red-shifted polarization signals ft chimoto et al.l [2007,) . The 
strength of the signal in the far wing, however, does not only 
depend on the velocity itself, but also on field inclination or gas 
temperature (Bellot Rubio 2009). 

In this contribution, we present measurements of line shifts 
in the wings of the Stokes I profile, which enables us to analyze 
the Doppler velocities, directly and quantitatively. 



2. Observation 

We used data o btained by the spectropolarimeter (SP) 
(iLites et al.l 1200 ll) of the solar optical telescope (SOT) 
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dTsuneta et alj l2008h onboard the Hinode satellite. The SP 
records the Stokes spectra of the two iron lines at 630.15 nm 
and 630.25 nm, Lande factors of g = 1.67 and g = 2.5, re- 
spectively. By scanning the spectrograph slit across the target 
in steps of 0."15 and achieving a wavelength sampling of 2.15 
pm/pixel, two dimensional maps of the field of view (FOV) are 
obtained. Since the width of the slit is equivalent to 0."16 and 
the pixel size along the slit is 0."16 on average, normal SP scans 
provide a spatial resolution of 0."32 dCenteno et al.l 120091) . For 
an exposure time of 4.8 s per slit position, the Icr noise level 
in the Stokes spectra is of the order of 10"^ Tc where Ic is the 
continuum intensity. 

Table 1. Data used in this study. From left to right: Sample name, 
NOAA active region, date of observation, heliocentric angle of 
penumbra, or field of view in quiet Sun samples. 



Data Sample 



Name 


NOAA 


Observation 






Spot A 


10923 


Nov 14"" 


2006 


0.982 - 


0.996 


SpotB 


10923 


Nov 14"" 


2006 


0.980 - 


0.994 


SpotC 


10930 


Dec U* 


2006 


0.991 - 


0.999 


SpotD 


10933 


Jan 05* 


2007 


0.996 - 


1.000 


QS 1 




Mar 10"' 


2007 


0.970 - 


1.000 


QS2 




Sep oe* 


2007 


0.988 - 


1.000 



For our study, we obtained the raw data summarized in 
Table [1] and reduced it, using the IDL routine 'sp_prep.pro' pro- 
vided with the 'Solar Soft' package. From the Stokes spectra, 
we computed maps of continuum intensity, I^, total polarization 
Ptot = / [(Q^ + + V2)/l2] i/2d^ and Doppler velocity along the 
LOS, Vdop. Furthermore, we increased the contrast of the maps of 
Ic and Ptot by applying a standard IDL high-pass filter to the pic- 
tures and then convolved the result with the original. As a conse- 
quence, the picture appears sharper, and it becomes easier to dis- 
tinguish between bright and dark features in the penumbra (e.g. 
the lateral brightening and the dark core of a filament). While 
the calculation of the maps of L and Ptot is straightforward, the 
derivation of Vdop requires some additional explanation. 

3. Method 

Since Hinode spectra lack telluric lines, it is not possible to use 
the latter to calibrate the velocity scale absolutely. To overcome 
this problem, we used an average QS profile for calibration pur- 
poses. The disadvantage of this method is, however, that the re- 
spective lines are displaced from laboratory wavelengt h by sev- 
eral e f fects, includ ing the convective blue-shift (CBS) dDravinsI 
119821) . iBeckl (|2005|) u sed the SIR code with a two-c omponent 
model atmosphere of iBorrero & Bellot Rubid (|2002|) to deter- 
mine the CBS at disk center. The corresponding velocities of the 
line core of the average QS profile of Fe 1 630. 15 and Fe 1 630.25 
were calculated to be -185 m s"' and -262 m s"' respectively 
(negative velocity corresponds to blue-shift). 

3.1. Velocity Calibration Using Quiet Sun Profiles 

We have argued that the line core of an average Stokes I profile 
of the QS is suitable for an absolute wavelength calibration, if 
it is corrected for the CBS. However, two effects, namely the 
moat flow and the influence of the magnetic field, have to be 
considered in the surroundings of sunspots. 



Due to the resid ual heliocentric angle, the radi ally outward 
oriented moat flow (Vargas Dommgu ez et al.ll2008h would lead 
to a systematic offset in the calibration, if only the Stokes I pro- 
files of one side of the spot vicinity were taken into account. 
Additionally, a magnetic field lifts the degeneracy of the atomic 
levels, which causes the profile to split and mak es it difficult t o 
determine the position of the line core correctly (Solanki! ! 19931) . 

Therefore, we decided to average all profiles of the QS sur- 
rounding the sunspot, but to exclude those profiles that exhibit 
a split. In order to meet this requirements, we created a mask 
which omits the entire spot and all pixels with significant mag- 
netic field strength. For these purposes, we used Ptot as a proxy 
and defined a cut-off criterion Ptot > 0. 1 • max(Ptot) for all pixels 
whose Stokes I spectra should not be averaged. As an example, 
the mask of Spot D is indicated in white in Fig.[T] 

Finally, the line minima of Fe I 630.15 and Fe I 630.25 in 
the average Stokes I profile of the QS are determined by fitting a 
2"*^ order polynomial to the core of the respective line (cf. dotted 
line in Fig.|2]i. 




Distance [arcsec] 

Fig. 1. Total polarization of data sample Spot D. The spectra of 
the pixels marked in white were not considered for an averaged 
Stokes I profile of the QS. The arrow points towards disk center 
at(x,y) ~ (87", 195"). 



3.2. Velocity Calibration Using Umbral Profiles 

iBeckersI (11977 1) showed that there is no vertical plasma motion 
in umbrae of sunspots at disk center Consequently, an average 
umbral profile should be red-shifted when compared to the line 
core of an average QS profile. To substantiate the va lues given by 
Beck (200i|), we followed a procedure described in lRezaei et al.l 
(2006). 

Due to the strength of the magnetic field, the core of umbral 
Stokes I profiles splits, and it is difficult to determine the core 
position correctly. Therefore, we determined the center-position 
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- midways between the two lobes - of umbral Stokes V profiles, 
which corresponds to the line core of the respective Stokes I pro- 
file. This parameter can be calculated with reasonable precision 
if molecular blends are not too strong and both lobes are well 
defined. If additionally LTE can be assumed, and no velocity 
gradients are present in the line forming region, e.g. above the 
umbra. Stokes V profiles are strictly antisymmetric with respect 
to the center-position (Auer & Heasley 1978). Thus, the center- 
position of antisymmetric Stokes V profiles defines a frame of 
rest on the solar surface. 

As a measure for the degree of antisymmetry, we used the 
amplitude asymmetry. It is defined as the difference of the abso- 
lute amplitudes of the two lobes like 



with ar and ab being the amplitude of the red and the blue lobe, 
respectively. We used a parabola fit around the extrema of both 
lobes to determine its amplitude. 

Following the arguments above, we considered only umbral 
Stokes V profiles with a \6a\ < 0.01 and determined their center 
position. The results are summarized in Table |2] 

Table 2. From left to right: Data sample, CBS of the average QS 
if the umbra is assumed to be at rest, and RMS values the center 
position of umbral Stokes V profiles with |^a| < 0.01. 




630.10 630.15 630.20 630.25 630.30 
Wavelength [nm] 



Fig. 2. Stokes I and Stokes V profiles (soUd line, measurements 
plotted with 'X' symbols) from a down-flow region (marked by 
the cross hairs in the inlet) in Spot A as well as the average 
Stokes I profile from the QS (dotted line, measurements indi- 
cated by symbols). The mean value of the bisector (solid) of 
the line wing (shaded gray) is used to derive the Doppler-shift of 
that pixel. Ic (dashed-dotted) is calculated by averaging the mea- 
surements below the two bold lines on the left and right sides of 
the spectrum. 



Wavelength Calibration 



Name 


CBS [m s-'] 


RMS Center Position [m s 




Fe I 630.15 Fe I 630.25 


Fe I 630.15 


Fe I 630.25 


Spot A 


-231 -296 


93 


83 


SpotB 


-245 -315 


94 


82 


SpotC 


-201 -295 


100 


92 


SpotD 


-231 -344 


147 


107 



'] 



We find that in all our data samples the average QS is in- 
deed blue-shifted with respect to the average center position of 
umbral Stokes V profiles obeying \6a\ < 0.01. Within the listed 
RMS errors (cf. Tab. |2]l, the values for the CBS agree with the 
ones reported by Beck (2005). Furthermo re, we can confirm t he 
applicabiUty of the procedure reported bv lRezaei et al.l (l2006l) . 

3.3. Flow velocity in the deep photosphere 

As we are interested in the plasma flow in the low photosphere, 
it is necessary to compute the line position as close to the con- 
tinuum as possible. On the one hand, the line forming region of 
Fe 1 630.15 encompasses higher atmospheric layers if compared 
to Fe I 630.25, but on the other hand, due to the lower Lande g 
factor, the former is less sensitive to disturbances by magnetic 
fields and instrumental crosstalk. In Fe I 630.15, a blend in the 
red wing close to the continuum obscures the profile and makes 
it impossible to derive a reliable bisector above 0.9 ■ Ic (cf. dot- 
ted line in Fig. |2]l. Similar obstacles occur in the spectra of Fe 
I 630.25 originating from the innermost penumbra. In this area, 
the magnetic field, and thus the line broadening, is so strong that 
the adjacent iron line serves as a blend, leaving the blue part 
of the continuum not defined above 0.9 ■ Ic. Therefore, we re- 
stricted the determination of the bisector of each line to a band 
of 0.7 < Ic < 0.9 (gray shaded region in Fig.|2]l. 



Next, we determined the bisector of both lines with sub-pixel 
accuracy by interpolating linearly between the spectral measure- 
ments. From that bisector (solid line in the middle of the shaded 
area in Fig. |2]l, we obtained an average spectral position of the 
line wing. Any offset between the latter position and the posi- 
tion of the line core of the averaged QS profile is interpreted as 
a Doppler-shift caused by the velocities of the plasma. Finally, 
the Doppler maps (cf. Fig.|5]l are obtained by computing the ve- 
locities that correspond to the Doppler-shifts at any given spatial 
position. 

As mentioned above, the low photospheric temperature of 
the umbra allows the formation of molecules. Their absorption 
bands mask the two iron lines to an extent where it becomes 
unreliable to compute a bisector Therefore, the Doppler velocity 
of pixels with an Ic of less that 33 % of that of the average QS 
were manually set to zero (cf. Fig.|5]l. 

3.4. Error discussion 

The biggest challenge in the calibration process is the determina- 
tion of the line core position of the average QS profile. Due to the 
finite spectral resolution, we have seven spectral measurements 
to fit the core. Using, for example, only five measurement to fit 
the core in sample 'Spot A' changes the CBS of the QS by -1-40 
m s In the QS significant magnetic fields are predominately 
found in the intergranular down-flow regions. As a result, we 
exclude more pixels showing down-flows and introduce a bias 
to our calibration (cf. Fig.[T]i. This bias leads to an overall blue- 
shift of the maps of the order of -10 m s"' to -30 m s"' in all 
data samples. The thermal drift in spectral direction is corrected 
by the 'sp_prep.pro' calibration routine. In all our data samples, 
it is within the range of +50 m s ' . Furthermore, we have shown 
that the values of the CBS of the QS inferred from model calcu- 
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lations and the values obtained from measurements differ by not 
more than +80 m s"'. 

Due to the difficulties mentioned above (especially the un- 
certainties in the position of the line core of the average QS) 
and since our fit procedure allows us to determine the line core 
within an accuracy of one tenth of a pixel, we assume a preci- 
sion of +0.1 km s"' for our wavelength calibration. As we do not 
consider any velocity signal -0.1 kms ' < vjop < 0.1 kms"' in 
the following, we are confident that any up- or down-flow is an 
actual velocity signal. 



down-flow shows peak velocities of up to 5.0 km s \ which are 
generally located at the penumbra-QS boundary, while weaker 
down-flows (Okms"' < Vdop < 2 kms"') appear at all radial dis- 
tances. The maximal up-flow velocities found in the other sam- 
ples are very similar. The maximum down-flow velocity, how- 
ever, shows variation of up to 80 % (cf. Table O. Interestingly, 
the former seems to change on a temporal scale of hours. In Spot 
A it peaks around 9.0 km s ', 9 hours later in Spot B it peaks at 
6.9 km s"'. Note that on a spatial average the penumbra shows 
a significant red-shift corresponding to the down-flow velocity 
listed in Table [3] 



4. Results 

4.1. Differences between Fe I 630. 15 and Fe I 630.25 

The resulting velocities derived from the Doppler-shifts of the 
two line wings are comparable: The RMS values of the maps 
deviate less than 0.17 km s ' . In general, the maps deduced from 
the wing of Fe I 630.15 show an overall blue-shift with respect 
to the maps deduced from the wing of Fe I 630.25. 

It seems that in all cases with significant differences 
(|vdop630.i5 - Vdop630.25l > 0.5 km s"'), the Fe I 630.25 line ex- 
hibits Doppler-shifts closer to the continuum than the Fe 1 630. 15 
line. This behavior is shown for a selected pixel with a large 
velocity discrepancy (Avjop = 1.5 km s ') in Fig. |2] Here the 
Stokes I and Stokes V profiles from a down-flow region in the 
outer penumbra of Spot A are depicted. The average Stokes I 
profile from the QS is indicated by the dotted line. 

We see a significant difference in the asymmetry of both 
lines, especially close to the continuum, which has a strong im- 
pact on our calculated velocity. In this example, the red wing 
of Fe I 630.15 is strongly shifted above 0.35 ■ Ic. The respec- 
tive wing of Fe I 630.25, however, shows a strong red-shift only 
above 0.77 • I^. Thus, the resulting velocities show significant 
differences. While the Doppler-shift of the average line wing 
(0.7 < Ic < 0.9) of Fe I 630.15 corresponds to 7.8 km s"', the 
respective Doppler-shift of Fe I 630.25 leads to 6.3 km s that 
is difference of 1.5 km s"'. In some cases, the asymmetries are 
above 0.9 ■ Ic and remain undetected in Fe I 630.25. Note, how- 
ever, that only the amplitude of the velocities, but neither the 
morphology of up- and down-flow patterns nor the direction of 
the flow are affected by this deviations. 

We do not have an explanation for these differences, but we 
are tempted to ascribe the variations in the maps of vjop to dif- 
ferences of the respective line parameters. Both lines have differ- 
ent g-factors, with Fe I 630.25 being more sensitive to magnetic 
fields. The line core of Fe I 630.15 forms at log t = -2.9, while 
the core of Fe I 630.25 forms around logr = -2.0 and there- 
fore encompasses lower atmospheric layers. As the strength of 
the magnetic field decreases with hight, its influence on the Fe 
I 630 .25 line is even more significant dKhomenko & ColladosI 
l2007h . 



Spot A 



Spot B 




0.5 1.0 1.5 
IDoppler velocityl [km/s] 



0.5 1.0 1.5 2.0 
IDoppler velocityl [km/s] 



SpotC 



Spot D 




0.5 1.0 1.5 
IDoppler velocityl [km/s] 

QS A 




0.5 Ifi^L 



0.5 1.0 1.5 
IDoppler velocityl [km/s] 

QSB 




0.5 1.0 1.5 2.0 
IDoppler velocityl [km/s] 



0.5 1.0 1.5 2.0 
IDoppler velocityl [km/s] 



Fig. 3. Histogram of binned Doppler velocities in sunspot (top 
rows) and QS (bottom row) samples in which the fraction of 
the penumbra or QS is depicted vs. the absolute value of the 
Doppler velocity. Black indicates up-flow, gray is down-flow. In 
all these measurements, the error is ±0.1 km s ' and indicated 
by the black error bars. 



4.2. Global Velocity Field 

As all data is taken almost at disk center, the maps of Doppler ve- 
locities (see Fig.|5]as an example) are dominated by plasma flows 
normal to the surface. We attribute any asymmetry between the 
center-side penumbra, which is overall blue-shifted, and the limb 
side, which shows an overall red-shift, to a combination of the 
horizontal EF and projection effects caused by the small helio- 
centric angle of the respective data sets. The up-flow in Spot D 
appears predominately but not exclusively in the inner penum- 
bra. It encompasses peak velocities around -2.0 km s"'. The 



The histograms in Figure [3] show the distribution of veloci- 
ties in the data sets in more detail. The absolute value of a ve- 
locity range of 10 m s ' is binned into one measurement and 
plotted on the abscissae. The fraction of the penumbra showing 
the respective velocity range is depicted on the ordinate. Down- 
flows appear in gray, while up-flows are drawn in black. Keep 
in mind that no unambiguous flow direction can be derived for 
|V(jop| < O.lkms"'. For reasons of comparison, QS data is evalu- 
ated in the same way (cf. bottom row of Fig. |3]l. 

In all sunspot data, the area occupied by up-flows is dom- 
inant at small velocities (|vdop| < 0.4 km s"'). The area show- 
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Fig. 4. Map of continuum intensity (Ic) of Spot D at 3°. White/black contours outline the inner/outer penumbral boundary and 
the arrow points towards disk center The three ellipses mark tracks in the penumbra at different radii for which Ic and vjop are 
investigated in detail in Sect. 14.31 The path is drawn counterclockwise and starts at the tip of the arrow. 



ing up-flows increases slightly and peaks at -0.1 < vjop < 
-0.2km s which seems to be a preferred velocity for up-flows. 
For larger velocities, this area decreases gradually and is equal in 
size to the down-flow area at velocities 0.4 < |vdop| < 0.6km s"', 
but negligible for Vdop < -1.5kms"'. The area showing down- 
flows, in contrast, decreases monotonously towards larger veloc- 
ities without any maximum and is the dominant one for |vdop| > 
0.6 km s"' . Even at |vdop| > 1 .5 km s"' , there is a significant area 
of the penumbra occupied with down-flows. 

In the QS data, the area occupied by up-flows is, however, 
always larger compared to the area showing down-flows, regard- 
less of velocity. Furthermore, the maximum down-flow velocity 
(3 km s"') is significantly smaller in the QS, while the maxi- 
mum up-flow velocity (-3 km s"' ) is about 50% larger as in the 
penumbra. 

These results are consistent with a fact mentioned above: 
Models and other observations imply an overall blue-shift of the 



QS, while our investigation yields a significant red-shift of the 
spatially averaged penumbra. 



4.3. Morphology of Velocity Field 

The Doppler map of Spot D shows a concentration of up-flows 
in the inner penumbra. They often appear as radially elongated 
patches next to each other separated by areas without any LOS 
motion. A typical up-flow area is about 2."4 long and 0."5 wide, 
resulting in a length/width ratio of about 5. Down-flows are pre- 
dominately visible at the penumbra-QS boundary. They are not 
elongated, but rather roundish. The length of the down-flow ar- 
eas ranges from l."3 to 3."0, while their width extends from 
l."3 to 2."0. Thus, with 1 to 1.5 the down-flow length/width ra- 
tio is smaller than for the up-flow areas. In order to investigate 
this morphology quantitatively and to compare the flow field to 
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30 

Distance [arcsec] 



Fig. 5. Same as Fig. HI but map of LOS velocities in Spot D as derived from the Doppler-shifts in the wing of Fe I 630.15 nm. 
Up-flows are shown from blue to green, down-flows are drawn from green to red. Note that the map is saturated at +1.5 km s and 
the velocities in the umbra are manually set to zero. 



the continuum intensity, we plotted Vdop and Ic along azimuthal 
paths. 

Fig. |6] shows Vdop and Ic in the inner penumbra. What is 
clearly visible are strong spikes of up-flow with velocities of up 
to -1.5 km s"' . Sometimes down-flow peaks are visible, but they 
are weaker and broader when compared to the up-flows. Areas 
of up-flow coincide well with areas of increased Ic, and areas of 
down-flow show a reduced Ic (cf. position B & A in the upper 
plot of Fig.|6]l. This is in agreement with the idea that up-flows 
transport hot plasma from lower layers to the surface, where it 
cools radiatively and appears bright. Due to the energy release it 
becomes more dense than the surrounding plasma and ultimately 
sinks back below the surface. On average, the plasma shows an 
up-flow of approx. -0.2 km s ' in the inner penumbra. 

In the mid penumbra (cf. Fig. |7]i, strong up-flow sites are 
seen less often and down-flows become more prominent. Here 
the penumbra still shows a net up-flow of about -80 m s"', and 



the correlation between up-flow and enhanced Ic is still valid 
(cf. F & G in Fig.|7j. Surprisingly, the strongest up-flow area of 
Spot D, that is -2.2 km s is located in the middle penumbra. 
With 2.3 km s ', the down-flows show velocities of comparable 
strength. 

In Fig. [8] Vdop and Ic of the outer penumbra are plotted. Here 
strong up-flow patches have almost disappeared and the domi- 
nant features are down-flows, with large velocity amplitudes of 
up to 3.2 km s"' . Note that not all up-flows are within the penum- 
bral boundary, but sometimes belong to granulation. Only one 
up-flow within the penumbra shows velocities of -1.6 km s"', 
while the rest is less strong. On average, the outer penumbra ex- 
hibits a down-flow of almost 0.45 km s"'. Interestingly, there is 
no distinct correlation between down-flows and a decrease of Ic 
in this part of the penumbra anymore. Position H & I in Fig.[8]in- 
dicate two strong down-flow places that show an Ic that is larger 
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Table 3. From left to right: Peak up- and down-flow veloci- 
ties derived from the Fe I 630.15 line wing, average velocity 
of penumbra in both lines wings, RMS values of differences be- 
tween the velocity derived from the two line wings. Negative 
values correspond to up-, positive values to down-flows. 



Mid Penumbra 



Doppler Velocities in Sunspot Penumbrae & Quiet Sun Data 
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Fig. 6. Continuum intensity (dotted) and Doppler velocity (solid) 
in the penumbra along the inner path shown in Fig. |4] and |5] 
Zero velocity is indicated by the solid horizontal line, while the 
dashed line marks the average velocity along that path. 



compared to the Ic of adjacent pixels which exhibit a weaker 
down-flow or even an up-flow. 

In the inner and mid penumbra, the flow pattern seems to 
follow the filamentary structure. This is not the case in the 
outer penumbra anymore, where the down-flow patches appear 
roundish and show no filamentary structure. In all azimuthal 
cuts, no distinct difference between the center and the limb side 
penumbra can be detected in the fl ow pattern, as it is the case in 
penumbrae away from disk center jTritschler et alj|2b04l) . 

4.4. Stokes V signal and magnetized flows 

The Stokes V profile from the down-flow region in, e.g. Fig. [2] 
shows strong asymmetries and more than two lobes. In order to 
explain the additional lobes with a simple scenario, we have to 
assume two magnetic components placed either above or next 
to each other in the resolution element. In both settings, one of 
these magnetic components has to be Doppler-shifted (in this ex- 
ample to the red) with respect to the other Thus, one component 
would represent a more or less steady magnetized background 
plasma, while the other component would incorporate strong and 
magnetized down-flows. 

The first scenario is more likely as it readily explains the 
shape of the Stokes I profile (cf. Fig.|2]l. The line wing is shifted 
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Fig. 7. 



Same as Fig.|6]but along the central path in Fig.|4] 

Outer Penumbra 




Fig. 8. Same as Fig. |6]but for the outer path in Fig. |4] 



towards the red, while the core is broadened but almost un- 
shifted. Therefore, it is probable that strong velocity gradients 
along the LOS are present. 

Keep in mind that in our case the velocities derived from the 
line wing of the two iron lines do underestimate the flow veloc- 
ity. Using the bisector to derive Doppler velocities will lead to 
average flow velocity in the resolution element. Assuming two 
components - one at rest, the other in motion - would automati- 
cally increase the velocities of the latter component. 
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4.5. Signatures of roll-type convection? 

In the framework of the gappy penumbral model energy is trans- 
ported by roll-type convection in the penumbral gaps, that is ar- 
eas without magnetic field below t = 1 . Thus, up-flows would 
be present in the center or the dark core of the filament, while a 
down-flow should be o bservable adja cent to it at the position of 
the lateral brightening (I Scharmer & S pruit 2006). Consequently, 
a down-up-down flow pattern ought to be visible on a scale of 
0.5" in cuts perpendicular to the filamentary structure of penum- 
brae. 

The plots of Fig.|6]to|8] however, show that strong up- and 
down-flow patches have a width of more than 0.5" and adjacent 
to them, the flow does not change its direction, but only eases in 
strength (cf. position B & D in Fig. |6]l. In rare cases, a reversal 
of the flow direction is observable on a scale of 0.5" (e.g. po- 
sition A & C in Fig. |6]l. The depicted flow pattern, however, is 
just opposite to the signature expected from roll-type convection. 
One of the very few examples, which shows a down-up-down 
flow pattern on a scale of less than 1" is indicated as E in Fig. 
|2](the peak-to-peak distance of the down-flow patches is 0.6"). 
A comparison of the maps of Ic and Vdop, however, shows that 
this flow pattern does not coincide with a (dark cored) filament. 
Other down-up-down flow signatures are observable, but have a 
width larger than 1" (e.g. the down-flow patches at position F 
show a peak-to-peak distance of 1.4"). 

A careful analysis of the other data samples does not yield a 
down-up-down flow pattern on scales of less than 1" that could 
be interpreted as a signature of roll-type convection in penum- 
bral filaments. 



5. Summary & Conclusions 

We presented two independent procedures to calibrate the wave- 
length scale of Hinode SP data. The first method uses the line 
core of an average Stokes I profile of the QS. This position was 
corrected for the CBS of the QS using the results teecfll2005h 
of model calculations (i Borrero & Bellot Rubioll2002l) . In a sec- 
ond approach, we assumed the umbra at rest and determined the 
center position, that is midways between the two lobes, of um- 
bral Stokes V profiles with an amplitude asymmetry of less than 
1%. This center position defines a frame of rest on the solar sur- 
face dRezaei et"al ]l2006h . Within the uncertainties, both methods 
yield the same results. 

In order to investigate the flow field in several sunspots close 
to disk center (0 < 9°), we determined Doppler velocities in the 
low photosphere by calculating an average bisector in the wings 
of both Fe lines. Even though in some cases the amplitudes of 
the flows derived from Fe I 630. 15 and Fe I 630.25 differ signif- 
icantly, the morphology of the flow field is not aff'ected by these 
differences. Furthermore, we find several prominent examples 
of multi-lobed Stokes V profiles which can be explained by a 
magnetized plasma in combination with strong velocity gradi- 
ence along the LOS. We therefore conclude that the downward 
component of the EE in the penumbra is magnetized. 

Since all the penumbrae we investigated are located close 
to disk center, we are able to visualize the LOS motion of the 
plasma normal to the solar surface, but cannot make a propo- 
sition conside ring the horizontal component of the EE. Other 
investigations dRezaei et al.ll2006 ?). however, show that at large 
heliocentric angle, the center position of Stokes V profiles is 
e.g. blue-shifted on the center-side penumbra, which indicates 
a magnetized horizontal component of the EE. 



Additionally, we find that the vertical flow pattern in sunspot 
penumbrae is structured on a small scale. Up-flows appear 
elongated and predominately in the inner penumbra. Down- 
flows concentrate at the penumbra-QS boundary and are rather 
roundish in shape. Ic and Vdop along azimuthal cuts in the penum- 
bra show that on average there is an up-flow present in the in- 
ner and mid-penumbra, while in the outer penumbra, on av- 
erage a down-flow is present. The up-flow patches in the in- 
ner penumbra and the down-flow locations in the outer penum- 
bra could be interpreted as the source s and the sinks of the EE 
dlchimoto & SOT/Hinode-teamll2009i) . even thoug h an attribu- 
tion t o individual flow channels is challenging (llchimoto et al.l 
^007h . 

Another finding of our investigation is the fact that the prop- 
erties of the penumbral flow field are strikingly different from 
that of the QS. In both iron lines, the area of the QS showing 
up-flows is always larger than the area showing down-flows, re- 
gardless of the velocity. This is not surprising as the up-flows are 
located in the granules and the down-flows are seen in the narrow 
lanes that form the intergranulum. Interestingly, in the penumbra 
up-flows dominate only at low velocities Ivjopl < 0.4 km s"'. At 
large velocities |vdop| > 0.6 km s"', the area occupied by down- 
flows outnumbers the one showing up-flows. This is not only in- 
dependent of the sunspot sample, but also irrespective of the Fe 
line we used to derive the Doppler velocity. Furthermore, in all 
penumbrae up-flows Vdop < -2 km s"' seem to be suppressed, 
while the velocity of up-flows in the QS can reach vjop = -3 
km s"' . On the contrary, down-flows with velocities of 9 km s"' 
are detected in some penumbrae; three times as large as in the 
down-flow regions of the QS. Additionally we find that on aver- 
age the penumbra shows a significant red-shift corresponding to 
a down-flow of more than 0.1 km s ' . 

If the penumbral energy transport is due to convective roll 
motion, a pronounced down-up-d own flow pattern should be ob- 
servable dScharmer & Spruiil2006.) in azimuthal cuts that cross 
(dark cored) filaments perpendicularly. In all our data sample, 
we do not find an indication for a roll-type motion of the plasma 
on a scale of the width of a filament. 

We can speculate why we do not see such a signature, and 
various explanation are possible. The easiest interpretation of 
this observation is that no roll-type convection exists in sunspot 
penumbrae. On the contrary, one could argue that we do not see 
hints of roll-type convection due to the limited spatial resolution 
of the Hinode SP. Another explanation could be that the Fe lines 
are not affected by the convective plasma motion because it is 
located only slightly above or entirely below t-\. 

The argument of limited spatial resolution always applies. 
Observations of (dark cored) penumbral filaments show that 
they have a width of less t han 1" so metimes with a dark core 
around 0.2" in wid th (Schar mer et al.ll2 002; Sutt erhn et al.l2004t 
iBellot Ru bio et al. 200 5J). In their MHD simulation of an entire 
sunspot, IRempel et al.i d2009t) found a down-up-down flow sig- 
nature with velocity variation of +1 km s"' across a filament. 
Even though the observation of the flow morphology in such 
dark cored penumbral filaments is at the limit of the spatial reso- 
lution of the Hinode SP, a distinct down-up-down flow pattern on 
scales of 0.5" should be observable in cuts perpendicular through 
the filament. If the plasma motion is entirely below r = 1, we 
have no chance to ever observe it. Yet, any plasma motion below 
T - \ will affect the higher atmospheric layers as well. So even if 
there is convective motion only below the surface, we should see 
at least a weak Doppler signal in the far wings of the Fe lines. 
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A more detailed investigation of the flow field in small scale 
features of the penumbra, e.g. bright penumbral grains, fila- 
ments, etc., will be addressed in a forthcoming publication us- 
ing the same data. Furthermore, it will be interesting to compare 
synthetic Stokes profiles computed from realistic simulations of 
radiative magneto-convection with our observations. 
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